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SECTION  1 


INTRODUCTION 


For  the  past  two  decades  optical  data  processing  (ODP)  has 
promised  a  vast  increase  in  processing  capacity  and  speed  over 
conventional  electronic  techniques.  This  promise  has  never  been 
fulfilled  for  several  reasons,  most  notably  because  of  the  lack 
of  a  practical  real-time  image  modulator,  or  light  valve,  and 
because  optical  techniques  were  almost  exclusively  limited  to 
linear  operations.  These  restrictions  have  been  removed  by  the 
development  of  the  liquid-crystal  light  valve  (LCVL)  by  Hughes 
Research  Laboratories  (HRL) ,  and  by  nonlinear  parallel¬ 
processing  techniques  developed  by  the  University  of  Southern 
California  (USC).  Thus,  it  is  important  to  determine  how 
successfully  non-linear  parallel-processing  techniques  can  be 
implemented  in  real  time  with  the  various  LCLVs. 

The  implementation  and  evaluation  of  these  techniques  have 
a  direct  relationship  to  current  Air  Force  technology.  Perti¬ 
nent  Air  Force  interests  include  multidimensional  real-time 
signal  and  image  processing  with  varied  applications,  including 
nonlinear  filtering  for  trajectory  control  and  guidance,  "smart" 
sensing,  picture  processing,  and  bandwidth  compression.  These 
technologies  could  benefit  substantially  from  the  increased 
processing  capacity  and  speed  that  the  proposed  research  may 
yield. 

Here  we  describe  nonlinear  optical  transformation  methods 
and  the  motivation  for  studying  a  real-time  application  of  those 
techniques.  In  Section  2  we  describe  the  progress  made  during 
the  current  program  year.  The  HRL-developed  photactivated 
twisted-nematic  LCLV,  the  variable  grating  mode  (VGM)  LCLV  modi¬ 
fication  and  its  system  configuration  for  optical  processing, 
and  optical  logical  operation  has  been  described  in  previous 
reports.  A  study  of  the  dynamic  performance  of  the  VGM  LCLV  has 
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been  performed  in  this  period  using  a  novel  technique  to  reduce 
the  time  of  VGM  domain  formation.  Empirical  modeling  of  the  VGM 
effect  has  also  been  performed. 

The  implementation  and  evaluation  of  these  techniques  has  a 
direct  relationship  to  Air  Force  interests  in  the  multidirec¬ 
tional  processing  of  real-time  signals  and  images,  with  applica¬ 
tions  including  nonlinear  filtering  for  guidance  and  trajectory 
control,  smart  sensing,  picture  processing,  and  bandwidth 
compression.  These  technologies  could  benefit  substantially 
from  the  increased  processing  capacity  and  speed  that  these 
techniques  may  yield. 

Until  now,  specified  nonlinear  operations  have  been 
performed  only  with  great  difficulty.  Coherent  optical  tech¬ 
niques  are  essentially  restricted  to  linear  opeations.  Digital 
processing  to  produce  nonlinear  transformations  is  possible,  but 
only  in  a  slow,  serial  fashion.  Certain  nonlinearities  can  be 
produced  by  special  photographic  techniques,  but  the  speed, 
accuracy,  reproducibility,  and  dynamic  range  of  these  techniques 
are  limited. 

We  have  been  pursuing  different  tasks  to  attempt  to  over¬ 
come  these  shortcomings.  The  first  made  use  of  special  half¬ 
tone  screens  to  modulate  the  input  image  in  conjunction  with 
coherent  optical  proessing  to  desample  in  the  Fourier  plane. 

This  technique  has  made  it  possible  to  implement  nonlinear 
effects  when  higher  orders  of  the  half-tone  diffraction  pattern 
are  examined  by  spatial  filtering.  Sawchuk  and  Dashiell  of  USC 
have  shown,  using  specially  fabricated  half-tone  screens,  how  a 
very  wide  class  of  two-dimensional  point  nonlinear  functions  can 
be  implemented  with  a  large  dynamic  range  as  a  function  of 
screen  design  and  diffraction  order.  The  nonlinearities  can  be 
continuous  or  discontinuous.  Operations  such  as  taking  loga¬ 
rithms,  exponentiation,  level  slicing,  intensity  bandstopping, 
and  histogram  equalization  can  be  performed.  We  have  expanded 
the  half-tone  screen  technique  by  substituting  a  real-time 


2 


photo-modulated  LCLV  for  the  static  photographic  recording 
medium,  and  we  have  successfully  demonstrated  a  logarithmic 
nonlinear  transformation  using  this  technique.  This  transforma¬ 
tion  is  useful  for  homomorphic-filtering  applications  as  we  have 
demonstrated.  In  cooperation  with  USC,  we  have  also  studied  the 
performance  potentials  and  limitations  of  this  implementation 
and  how  to  iteratively  modify  and  improve  the  LCLV  and  half-tone 
masks . 

A  second  general  method  which  we  are  presently  studying 
also  overcomes  the  limitations  of  serial  or  photographic  proces¬ 
sing;  it  employs,  in  one  realization,  a  liquid-crystal  effect  - 
VGM  -  that  can,  when  incorporated  into  a  new  type  of  photocon- 
ductive  structure,  automatically  map  image  intensity  variations 
into  positions  in  Fourier  space.  Filtering  and  reconstructing 
can  then  yield  many  desired  nonlinear  transformations  of  the 
image  without  the  need  for  specially  constructed  half-tone 
masks.  A  new  additional  parameter,  the  intensity,  has  thus  been 
made  available  for  optical  image  and  data  processing.  Recogniz¬ 
ing  that  logical  operations  are  merely  a  special  case  of 
non  L . near  operations,  we  have  demonstrated  a  unique  and  highly 
advantageous  optical  computing  scheme  using  the  VGM  technique. 
The  VGM  device  is  still  in  an  early  stage  of  development  and 
much  material  research  and  device  development  is  still  needed  to 
make  it  into  a  practical,  real-time,  reliable  optical  image 
modulator.  This  work  is  still  continuing.  Alternative  practi¬ 
cal  realizations  of  intensity  to  positional  mapping  are  also 
being  studied  that  may  overcome  some  of  the  shortcomings  of  VGM, 
yet  retain  the  great  flexibility  offered  by  intensity-to- 
positional  coding. 


SECTION  2 

PROGRESS  DURING  CURRENT  YEAR 

During  this  program  year,  we  pursued  several  tasks  which 
aimed  at  further  improving  our  understanding  of  the  VGM  effect 
and  at  developing  a  more  practical  VGM  photoact ivated  light 
valve  device.  In  the  VGM  effect,  certain  liquid-crystal  ele¬ 
ments  exhibit  a  phase  grating  whose  local  period  is  proportional 
to  the  locally  applied  voltage.  When  incorporated  into  a  light 
valve  structure  along  with  a  photoconductive  layer,  local  varia¬ 
tions  in  the  spatial  distribution  of  light  signal  or  image 
intensity  are  converted  into  local  variations  in  grating  period. 
In  the  Fourier  transform  plane,  a  new  parameter  -  the  intensity 
is  coded  into  the  spatial  distribution  in  a  way  that  allows 
great  freedom  to  perform  nonlinear  optical  data  and  image  pro¬ 
cessing,  and  to  perform  optical  computations  as  well.  In  this 
period,  we  continued  to  examine  many  aspects  of  the  device, 
including  the  temporal  dynamic  properties  of  the  light  valve 
using  novel  schemes  to  increase  the  speed  of  response.  The 
detailed  optical  polarization  properties  of  the  VGM  diffraction 
patterns  were  also  studied  with  the  goal  of  completing  the 
modeling  of  the  molecular  configuration  of  VGM.  We  also  studied 
alternative  schemes  to  the  VGM  wherein  the  main  feature  of 
converting  intensity  variations  to  positional  variations  would 
be  included  along  with  the  desired  faster  time  response 
characteristics . 

A.  IMPROVED  VGM  RESPONSE  TIME 


At  the  end  of  the  last  period  we  began  a  study  which  showed 
promise  of  improving  the  VGM  response  time  by  spoiling  the  long 
range  order  of  the  periodic  domains.  In  this  sub-section  we 
will  review  these  results  and  discuss  our  progress  during  this 
year . 


From  the  device  standpoint,  the  VGM  response  times,  as  we 
have  noted  before,  are  too  slow  for  many  applications.  It  would 
be  desirable  to  have  responses  compatible  with  TV  frame  rates  or 
faster,  on  the  order  of  tens  of  milliseconds.  In  this  section 
we  describe  a  novel  approach  to,  and  some  preliminary  results 
for,  achieving  faster  dynamic  response.  It  is  based  on  the 
intuitive  notion  that  the  VGM  domains  evolve  slowly  in  time, 
both  along  the  domain  direction  and  in  the  perpendicular 
direction. 

There  is  some  visual  evidence  for  the  slow  evolution  of  the 
domains.  When  the  voltage  is  changed  in  a  VGM  cell,  say  to  a 
slightly  higher  value,  under  the  polarizing  microscope  one  can 
see  the  new  additional  lines  coming  into  the  field  of  view.  New 
lines  form  by  "unzipping"  one  or  more  of  the  numerous  fork¬ 
shaped  dislocations  which  occur  throughout  the  area  of  the  VGM 
structure.  The  unzipping  process  typically  traverses  -20  m  at  a 
velocity  on  the  order  of  magnitude  of  10  m/sec. 

It  seemed  possible  that  if  the  length  of  the  VGM  grating 
domain  could  be  shortened  radically  so  that  the  unzipping  would 
not  have  to  travel  far,  that  the  time  of  formation  of  the  new 
VGM  steady  state  could  be  shortened.  One  way  to  cover  a  two- 
dimensional  area  with  short  VGM  domains  would  be  to  have  the 
longer  domains  interrupted  by  a  fixed  set  of  parallel  interrup¬ 
tions  perpendicular  to  the  length  of  the  domains.  The  interrup¬ 
tion  would  have  to  present  a  strong  electrical  and/or  mechanical 
perturbation. 

We  constructed  a  special  cell  with  one  narrow  1-mm  groove 
parallel  and  another  perpendicular  to  the  VGM  domains.  No 
difference  in  response  times  was  noted  between  these  grooves, 
nor  between  the  grooves  and  the  normal  large  area  VGM  cell 
geometry. 

Next  we  fabricated  a  set  of  holographically  produced  grat¬ 
ings  with  a  one  thousand  times  smaller,  1  m  period.  These 
gratings  were  formed  on  ITO  coated  glass  substrates  etched 
through  the  conductive  coatings.  The  conductive  strips  were  all 


electrically  connected  at  one  edge.  The  design  of  the  substrate 
is  shown  in  Figure  1 .  Two  orthogonal  gratings  were  put  on  each 
substrate  so  that  both  parallel  and  perpendicular  opposing  elec¬ 
trode  gratings  could  be  configured  in  the  complete  liquid  crys¬ 
tal  sandwich  cell  sturcture.  This  assured  that  both  these 
arrangements,  and  the  grating-free  areas  as  well,  could  all  be 
studied  under  the  same  conditions  of  thickness,  temperature 
impurities,  and  aging. 

The  VGM  domains  in  a  working  cell  6  m  thick  with  20  V-dc 
applied  are  shown  in  Figure  2.  The  domains  are  perpendicular  to 
the  two  1-  m  holographic  gratings.  In  the  photograph,  that 
grating  cannot  be  seen,  and  only  large  occasional  defects  along 
the  direction  of  the  grooves  show  up.  What  is  noteworthy  in 
this  typical  photograph  is  the  large  number  of  forked  disloca¬ 
tions,  as  compared  with  a  cell  made  of  normal  substrates  without 
any  gratings.  The  domains,  however,  are  not  merely  1  m  long 
(the  period  of  the  grating),  but  typically  extend  to  10  m  in 
length.  In  a  typical  cell  without  gratings  the  domains  may 
extend  several  millimeters  before  forking. 

The  LC  mixture,  Merck  NPV,  was  aligned  perpendicular  to  the 
direction  of  the  grooves  of  the  holographic  grating  in  those 
experiments  that  we  will  report  on  here.  In  another  set  of 
experiments  with  the  alignment  parallel  to  the  grooves,  no 
significant  effects  were  noted. 

The  temporal  response  characteristics  proved  to  be  complex 
and  impossible  to  describe  merely  with  a  single  set  of  charac¬ 
teristic  rise  and  decay  times.  The  responses  are  non-exponen¬ 
tial  and  multiprocess  with  a  pronounced  memory  or  hysteresis. 
There  is  also  a  secular  degradation  in  response.  All  of  these 
complexities  remain  to  be  analyzed. 

The  apparatus  shown  in  Figure  3  was  employed  to  perform  the 
response  time  measurements.  The  measurements  were  taken  on  the 
first  order  diffracted  light  from  the  VGM  cell.  The  diffracted 
spot  of  the  laser  beam  was  from  3  to  5  mm  in  diameter  at  the 
detector  20  cm  from  the  LC  plane.  The  detector  was  positioned 
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in  each  experiment  at  the  angle  corresponding  to  the  steady 
state  diffraction  angle  for  the  particular  field  applied.  A 
mask  was  positioned  in  front  of  the  detector  to  give  an  angular 
resolution  of  1  mm  at  20  cm  from  the  VGM  cell  plane. 

To  measure  the  temporal  response  of  this  cell,  the  voltage 
on  the  cell  was  repetitively  switched  from  one  dc  value  to 
another  in  a  rectangular  wave  pattern.  One  detector  was  placed 
at  each  of  the  diffraction  positions  appropriate  for  each  vol¬ 
tage.  The  outputs  of  these  detectors  were  displayed  on  dual 
traces  of  an  oscilloscope. 

For  example,  in  Figure  4  the  top  trace  records  the  output 
of  the  detector  placed  at  the  30-V  diffraction  position,  while 
the  bottom  trace  simultaneously  records  the  output  of  the  detec¬ 
tor  at  the  20-V  diffraction  position  as  the  voltage  is  repeti¬ 
tively  switched  from  20  to  30  V  in  a  2-sec  period.  The  motion 
of  the  diffracted  light  beam  from  one  detector  to  the  other  is 
complex.  Not  only  does  the  beam  travel  from  one  detector  to  the 
other,  but  it  also  undergoes  a  reformation  and  condensation  as 
it  approaches  steady  state,  making  the  temporal  response  very 
complicated,  as  can  be  seen  from  the  oscilloscope  traces.  There 
is  also  a  sudden  flash  of  light,  largely  centered  around  the 
current  position,  whenever  the  voltage  is  switched.  This  flash 
shows  up  as  a  strong  spike  on  the  traces.  From  Figure  4  we 
estimate  the  rise  time  to  be  600  t  200  msec  and  the  delay  time 
is  now  only  in  the  tens  of  milliseconds!  In  Figure  5  the  cell 
was  refilled  and  traces  were  recorded  from  30  to  40  V.  The  top 
trace  corresponds  to  40  V.  Here  the  rise  times  were  on  the 
order  of  100  msec,  but  the  decay  times  were  600  and  1,700  msec, 
respectively,  for  30  and  40  V.  This  was  clearly  not  a  reproduc¬ 
tion  of  the  remarkably  rapid  decay  observed  in  Figure  4.  At 
this  stage  in  our  experiments,  we  have  been  able  to  reproduce 
the  fast  performance  only  from  time  to  time  and  in  only  certain 
of  the  cells.  However,  we  are  greatly  encouraged  to  have  seen 
such  significant  improvement  in  temporal  response  using  this  new 
technique.  For  example,  after  the  gradual  deterioration  in  VGM, 


Figure  4.  VGM  temporal  response;  +20  to  +  30  V  dc 
1  urn  grating  substrates,  500  msec/cm. 


Figure  5.  VGM  temporal  response;  +30  to  +40  V  dc, 
1  um  grating  substrates,  2  sec/cm. 


performance  of  the  particular  cell  studied  in  Figures  4  and  5 
(noted  by  higher  thresholds  and  increasing  conductivity)  could 
no  longer  be  tolerated  and  the  cell  was  refilled.  The  decay 
times  were  four  times  slower  than  the  original  values.  Upon  a 
later  second  refilling,  the  times  were  even  slower.  This 
phenomenon  is  not  yet  understood. 

Next,  quasi-ac  experiments  were  tried.  If  operational,  ac 
would  avoid  electrochemical  poisoning  of  the  liquid  crystal. 
Switching  from  the  same  positive  to  negative  voltage  values, 
e.g.,  +20  V  to  -20  V  in  square  wave  fashion  also  produced  fast 
decay  times  of  the  order  of  50  msec  and  rise  times  of  the  order 
of  a  second.  But  as  the  diffraction  angle  is  the  same  for  the 
positive  as  for  the  negative  voltage,  and  since  there  is  some 
memory  of  previous  states  at  the  switching  times  we  emplopyed, 
we  abandoned  that  choice  of  equal  voltage  as  being  misleading. 

To  take  advantage  of  quasi-ac  and  yet  have  more  meaningful 
temporal  response  measurements,  we  repetitively  performed 
switchings  from  +30  to  -40  V.  Figure  6  illustrates  the  VGM 
response  of  a  1  :.m  holographic  grating  cell.  It  should  be  noted 
that  unless  otherwise  stated  the  configuration  normally  employed 
is  with  both  electrode  and  counter  electrode  gratings  parallel 

to  each  other.  In  the  top  trace  the  rise  time  from  +30  to  -40  V 

is  300  msec  and  the  decay  time  >10'  msec.  In  the  lower  trace, 
measured  on  the  second  detector,  the  rise  time  from  -40  to  +30 

is  1  sec  and  the  decay  time  from  +30  to  -40  is  *50  msec.  The 

asymmetry  in  rise  times  is  not  understood. 

The  production  of  holographic  gratings  is  time  consuming 
and  subject  to  many  difficulties  and  failure  modes  in  the  expo¬ 
sure  and  etching  phases.  It  seemed  that  it  might  be  possible  to 
spoil  the  long  range  order  in  the  VGM  domains  by  a  simpler 
method:  by  producing  linear,  randomly  spaced  scratches  on  the 

substrate  surface. 

One  method  of  liquid  crystal  alignment,  often  used  for  immedi 
ate  results,  is  rubbing.  The  substrate  is  rubbed  with  coarse  mate 
rial,  perhaps  diamond  paste,  to  effect  a  quasi-linear  pattern  of 
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VGM  temporal  response;  +30  to  -40V  dc 
1  ym  grating  substrates,  300  msec/div 


surface  forces  which  align  the  liquid  crystal.  There  is  no 
definite  period  in  these  linear  scratches.  They  are  randomly 
spaced.  We  first  rubbed  the  surface  in  one  direction  and  then 
cross-rubbed  it  in  the  orthogonal  direction.  The  liquid  crystal 
successfully  aligned  in  the  last  direction  that  was  rubbed,  but 
the  underlying  cross-rubbing  seemed  to  achieve  the  desired 
effect  of  breaking  up  the  long  range  VGM  domain  order.  This  can 
be  seen  in  Figure  7,  a  polarizing  microscope  photograph  of  a 
cross-rubbed  PVA  substrate  cell  spaced  at  6  m«*  and  filled  with 
Merck  NPV  with  50  V  dc  applied.  The  long  range  order  is  clearly 
spoiled.  The  temporal  behavior  of  this  cell  is  shown  in  Figure  8. 
The  top  trace  -40  to  +60  V  shows  a  rise  time  of  1.9  sec  and 
decay  time  of  =■  1  0  1  msec.  Even  this  crude  technique  has  signifi¬ 
cantly  reduced  the  response  time  of  the  VGM  device! 

A  set  of  0.4  pm  holographic  gratings  was  fabricated  and 
tested.  This  period  is  significantly  smaller  than  the  typical 
VGM  instability  period  and  it  was  unclear  whether,  for  best 
results  in  improving  the  response  time,  the  grating  period 
should  be  larger  or  smaller  than  the  VGM  period.  A  scanning 
electron  microscope  photograph  of  the  0.4-um  grating  of  ITO 
onglass  at  40,000  diameters  is  shown  in  Figure  9.  The  results 
with  this  grating  have  not  been  as  dramatic  as  with  the  1  um 
grating,  but  a  significant  improvement  over  no  grating  at  all 
was  still  observed.  In  Figure  10  we  display  the  results  for 
perpendicularly  disposed  grating  electrodes  and  for  +40  to 
+30  V  dc.  The  rise  and  decay  times  are  ~400  msec.  This  was 
typical  of  the  usual  grating  arrangement  as  well. 

A  6-um  set  of  holographic  grating  substrates  were  fabrica¬ 
ted  into  VGM  cells  and  evaluated.  This  6-um  period  is  somewhat 
below  the  typical  period  of  the  normal  VGM  phenomena.  With 
these  substrates  we  could  not  produce  clear  VGM  phenomena.  The 
diffraction  was  extremely  smeared  and  confused. 

To  illustrate  the  frequency  response  of  the  VGM  phenomena 
in  the  1  pm  grating  cells,  we  show  in  Figure  11  (a),  (b),  and 
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Figure  7.  VGM  domains  in  cross-rubbed  PVA  cell,  180 
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Figure  8.  VGM  temporal  response;  -40  to  +60  V  dc, 
cross-rubbed  PVA  substrates.  2  sec/div 


12778-4 


Figure  9.  SEM  photograph  of  0.4  urn  holographic 
grating  of  ITO  on  glass. 
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Figure  10.  VGM  temporal  response;  +40  to  +30  V  dc 
0.4  um  grating  substrates  perpendicu¬ 
larly  disposed,  2  sec/div. 
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1.  VGM  frequency  response;  -40  to  +20  V, 

1  um  grating  substrates,  1  Hz  to  30  Hz 
sweep.  (a)  Square  wave;  (b)  Sine  wave 
(c)  Sawtooth  wave. 


(c),  the  -40  to  +30  V  temporal  response  using  square,  sine  and 
sawtooth  wave  excitation,  respectively,  swept  from  1  to  30  Hz. 
There  is  clearly  significant  response  to  beyond  10  Hz.  Most 
noteworthy  is  the  strong  spiking,  mentioned  earlier,  which  can 
now  be  seen  to  be  most  prominent  at  a  certain  intermediate  range 
of  driving  frequencies,  depending  upon  whether  square,  sine,  or 
sawtooth  waves  are  employed.  It  would  seem  possible  to  drive 
devices  with  these  waveforms  and  implement  the  rapid  response 
implied  by  these  experiments. 

B.  PHYSICAL  ORIGIN  OF  THE  VARIABLE  GRATING  MODE  EFFECT 

Obviously,  the  important  VGM  device  operational  properties 
depend  on  the  detailed  nature  of  the  grating  formation  and 
dynamic  orientation  process.  We  have,  in  conjunction  with  our 
colleagues  at  USC,  made  some  experimental  and  theortical  studies 
this  year  to  elucidate  the  fundamental  nature  of  the  grating  and 
of  the  physical  mechanisms  that  give  rise  to  the  observed 
periodic  instability  and  its  dynamics. 

The  grating  produced  by  the  periodic  spatial  re-orientation 
of  the  nematic  liquid  crystal  molecules  is  quite  unusual  and 
gives  rise  to  striking  polarization-dependent  properties  in 
optical  diffraction.  By  studying  the  intensity  and  polariza¬ 
tion  of  the  diffracted  light  for  each  diffracted  order,  we  have 
been  able  to  determine  the  spatial  distribution  of  the  molecular 
orientation  within  the  VGM  liquid  crystal  layer.  The  polariza¬ 
tion  dependence  of  the  diffraction  is  due  to  the  birefringent 
phase  grating  in  which  the  orientation,  the  principal  axis  of 
the  index  ellipsoid,  varies  periodically  both  in  the  plane  of 
the  grating  and  normal  to  the  plane.  The  spatial  distribution 
of  the  ends  of  the  liquid  crystal  molecules  is  approximately 
cycloidal.  This  was  determined  by  matching  the  experimental 
results  with  theoretically  calculated  values  for  different 
possible  models.  The  technical  details  are  described  in  a 


L9 


preliminary  publication  reproduced  in  the  appendix.  A  more 
detailed  publication  is  being  prepared  as  well. 

A  theoretical  study  of  the  spatial  distribution  of  the 
liquid  crystal  molecules  is  also  being  carried  out  by  direct 
minimization  of  the  free  energy.  The  preliminary  results  of 
this  study  show  that  the  angular  dependence  observed  experimen¬ 
tally  can  be  predicted  theoretically.  The  origin  of  the  VGM 
effect  is  still  in  doubt,  but  the  evidence  strongly  points  to 
the  converse  f lexoelectric  effect.  The  details  of  this  theoret¬ 
ical  work  to  date  are  also  reproduced  in  the  appendix. 

The  dynamics  of  grating  formation,  reorientation  and  relax¬ 
ation  are  subjects  of  on-going  experimental  and  theoretical 
investigation.  Understanding  of  the  basic  physical  principles 
underlying  these  effects  is  vital  to  the  success  of  efforts  to 
improve  the  VGM  response  time. 

C.  ALTERNATIVES  TO  THE  VGM  DEVICE 

The  variable  grating  mode  liquid  crystal  effect  provides 
only  one  possible  means  of  achieving  parallel  intensity  to  posi¬ 
tion  encoding.  The  processing  potential  of  the  intensity-to- 
position  algorithm  is  so  great  that  we  have  been  exploring  alter¬ 
native  implementations.  One  promising  possibility  is  electro¬ 
optic  beam  deflection,  as  discussed  previously.  The  angle  that 
light  is  deflected  by  refraction  from  a  prism  of  a  small  angle  is 
approximately  proportional  to  the  prism  angle  and  the  index  of 
refraction.  The  index  refraction  can  be  varied  y  electro- 
optically  in  many  materials;  the  effect  is  particularly  large  in 
LCs  where  the  change  in  index  between  the  extremes  of  ordinary  and 
extra-traordinary  polarizations  can  be  as  much  as  0.2  with  rela¬ 
tively  small  voltages  (~ 1 0  V).  By  varying  the  voltage  locally  to 
conform  to  a  2-D  image  which  illuminates  a  photoconductive  layer  on 
top  of  a  prismatic  wedge  of  LC,  one  would  have  a  device  that 
transforms  local  image  intensity  into  local  angles  of  deflection. 


In  the  back  focal  plane  of  a  lens,  these  angles  are  (napped  into 
displacements,  i.e.,  a  spatial  Fourier  transform  space,  with  the 
new  additional  parameter  of  intensity  included  as  well,  just  as 
in  the  VGM  device.  In  the  last  period  we  demonstrated  electro¬ 
optic  deflection  with  liquid  crystal  material.  The  problem 
encountered  was  the  very  slow  response  time  of  thick  liquid 
cells.  In  order  to  bring  the  response  time  to  the  millisecond 
region  we  conceived  of  an  array  of  replicated  little  deflecting 
prisms.  Since  only  the  wedge  angle  determines  the  deflection  of 
the  prism  and  not  its  thickness,  the  array  of  little  prisms 
could  effect  the  same  deflection  without  the  disadvantages  of 
slow  speed  and  cloudiness  caused  by  large  thicknesses.  Our  past 
attempts  to  construct  such  structures  in  metal  and  in  glass  have 
thus  far  proven  difficult  and  expensive.  However,  we  have 
devised  a  new  idea  to  solve  the  problem.  We  are  in  the  process 
of  fabricating  the  periodic  structure  by  using  a  technique 
developed  for  other  purposes  at  our  laboratory:  the  anisotropic 
etching  of  silicon  for  the  production  of  blazed  deffraction 
gratings.  This  work  was  done  to  much  finer  tolerances  and 
smaller  dimensions  than  we  require  for  our  application. 

Anisotropic  etching  of  single  crystal  silicon  is  crystal 
orientation  dependent.  Silicon  has  a  diamond  cubic  crystal 
structure.  The  anisotropic  etch  rate  in  the  (100)  direction  is 
typically  greater  than  the  (110),  which  is  greater  than  the 
(111)  direction.  The  etching  process  is  influenced  by  the 
crystallographic  orientation  of  the  substrate,  the  alignment  of 
the  grating  mask  pattern  to  the  crystal  orientation,  and  the 
etchant  chemical. 

The  fabrication  process  basically  involves  a  substrate 
with  the  proper  crystal  orientation  to  achieve  a  particular 
grating  shape  or  blazed  angle.  A  Si02  layer  is  thermally  grown 
onto  the  substrate  to  serve  as  a  chemical  mask  for  the  aniso¬ 
tropic  etching.  Lithographic  patterning  of  the  grating  is 
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initially  done  in  a  photoresist  layer  over  the  oxide.  Alignment 
of  the  grating  pattern  to  the  crystal  orientation  is  important 
to  achieve  proper  groove  shape.  This  pattern  is  transferred 
into  the  oxide  coatig  by  acid  etching  or  by  reactive  plasma 
etching.  The  substrate  is  then  exposed  to  the  anisotropic 
chemical  etchant,  such  as  potassium  hydroxide  (KOH). 

Although  symmetric  grating  profiles  can  be  fabricated  using 
this  technique,  well  controlled  blazed  profiles  are  also 
achievable.  For  example,  if  a  Si02  mask  is  aligned  along  a 
(110)  direction  of  a  (100)  silicon  slice  (see  Figure  12),  V- 
shaped  grooves  are  produced  with  smooth  (111)  walls.  Unlike 
isotropic  chemical  etching,  the  etch  progresses  in  the  (100) 
direction  and  is  terminated  at  a  depth  of  .707  of  the  oxide 
opening  width,  where  the  etch  front  meets  the  (111)  planes, 
intersecting  the  (100)  plane  at  the  edge  of  the  mask  opening. 

The  apex  angle  of  the  groove  is  109.47".  To  fabricate  blazed 
structures,  the  silicon  substrate  can  be  cut  and  polished  to 
achieve  the  desired  blazed  angle  (see  Figure  13).  Note  that 
apex  angles  of  either  109.47’  or  70.53*  are  attainable, 
depending  upon  the  crystal  cut.  Recently,  we  have,  on  another 
program,  fully  fabricated  a  short  period  (two  micrometers) 
grating  in  a  single  crystal  silicon  with  significantly  reduced 
flat  top  defect  (Figure  14).  This  period  is  several  times 
smaller  than  that  which  would  be  required  for  this  program. 
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I.  INTRODUCTION 

A  wide  variety  of  one-  and  two-dimensional  operations  are  necessary 
tor  lull-scale  implementation  of  parallel  optical  processing  and 
computing  systems.  Incoherent-to-coherent  conversions  are  often 
required  for  algorithms  involving  spectrum  analysis  and  modifi¬ 
cation.  correlation,  convolution,  and  holographic  image  formation, 
particularly  when  the  information  to  be  processed  is  available  only  in 
time-sequential  or  matrix-addressed  raster  format.  A  number  of 
one-  and  two-dimensional  spatial  light  modulators  capable  of  this 
type  ot  image  transduction  are  described  within  this  special  issue  of 
Optical  Engineering,'-*  as  well  as  in  several  review  articles. 1_i 

Other,  equally  important  processing  and  computing  functions, 
such  as  logic  operations,  programmable  matrix  addressing,  binary 
addition,  linearity  compensation,  and  input-output  nonlinearities 
le  g.,  exponentials,  logarithms,  power  laws,  thresholds,  level  slices, 
and  level  restoration),  have  proven  particularly  difficult  to  implement. 
All  of  these  functions,  on  the  other  hand,  can  be  implemented  by 
means  of  some  form  of  intensity-to-position  encoding  in  conjunction 
with  either  fixed  (single  function)  or  programmable  (multifunction) 
masks  This  general  statement  follows  from  the  realization  that  all  of 
the  functions  listed  above  are  special  cases  of  data-dependent  multi¬ 
plications.  in  which  the  input  value  (e  g.,  pixel  intensity)  selects  the 
appropriate  multiplier  (e  g.,  mask  location)  to  obtain  the  desired 
product  (e  g.,  output  intensity). 
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The  variable  grating  mode  liquid  crystal  device  (VGM  LCD)’"12 
transforms  input  intensities  to  spatial  positions  when  used  in  conjunc¬ 
tion  with  a  Fourier  transform  lens.  The  nature  of  this  image  trans¬ 
formation  can  be  realized  in  the  following  manner.  The  VGM  LCD 
primarily  consists  of  a  photoconductive  layer  in  series  with  a  layer  of 
nematic  liquid  crystal  mixture.  A  dc  bias  voltage  is  applied  across  the 
device  to  provide  a  voltage  division  between  the  two  layers.  Within  a 
given  image  pixel,  the  input  intensity  decays  the  voltage  across  the 
photoconductive  layer  and  correspondingly  enhances  the  voltage 
across  the  liquid  crystal  layer.  The  photoconductor  thus  implements 
an  intensity-to-voltage  conversion.  The  nematic  liquid  crystal  mix¬ 
ture  employed  in  the  device  has  the  unusual  property  that  the 
alignment  of  the  liquid  crystal  molecules,  which  is  homogeneous  in 
the  quiescent  state,  exhibits  spatially  periodic  modulation  when  a 
bias  voltage  is  applied  across  the  layer.  This  modulation  results  in  a 
birefringent  phase  grating1 5  characterized  by  a  spatial  frequency  that 
depends  linearly  on  the  applied  voltage.  The  effect  of  the  liquid 
crystal  layer  is  thus  to  implement  a  voltage-to-spatial  frequency 
conversion.  If  both  layers  are  considered  together,  the  entire  device  is 
thus  seen  to  perform  an  image-wise  intensity-to-spatial  frequency 
conversion,  which  can  be  modified  to  the  more  general  intensity-to- 
position  transformation  by  placing  a  Fourier  transform  lens  behind 
the  VGM  LCD.  Collimated  readout  illumination  normally  incident 
on  the  device  (at  a  wavelength  of  photoconductive  insensitivity)  is 
angle  encoded  within  each  image  pixel  by  diffraction  from  each 
induced  phase  grating  and  subsequently  angle-to-position  mapped 
by  the  Fourier  transform  lens  into  its  focal  plane. 

This  type  of  process  is  shown  schematically  in  Fig.  1 ,  in  which  the 
input  image  is  assumed  to  consist  of  two  separate  regions  of  differing 
intensity.  The  VGM  LCD  encodes  both  regions  with  different  spatial 
frequencies,  resulting  in  separated  diffraction  orders  in  the  filter 
(Fourier)  plane  Insertion  of  an  appropriate  spatial  filter  or  pro¬ 
grammable  mask  (not  shown  in  Fig.  1 )  into  the  Fourier  plane  allows 
the  separated  orders  to  be  selectively  modified  to  implement  any 
desired  data-dependent  multiplication  or  point  nonlinearity.  In  the 
reconstructed  output  image,  all  regions  of  equal  input  intensity  are 
modified  identically,  irrespective  of  their  location  in  the  input  image 
field.  Thus,  all  of  the  data-dependent  multiplications  are  performed 
in  parallel.  Functional  programmability  is  achieved  by  replacement 
or  reprogramming  of  the  Fourier  plane  mask,  which  need  only  be  a 
low  resolution  device  with  a  total  number  of  resolution  elements 
equal  to  the  number  of  gray  levels  required  to  be  processed. 

The  overall  input-output  characteristic  of  nonlinear  function 
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«ion  in  a  VGM  LCD. 


Fig.  2.  VGM  noniinaar  procaaaing.  Tha  ovarali  input-output  charact aria- 
tic  can  bo  found  by  (tapping  through  tha  auccaaaivo  noniinaar  tranafor- 
mationa.  including  (1)  tha  intanaity-to-apatial  fraquancy  convaraion.  (2) 
apatial  fiitaring.  and  (3)  intensity  d  at  act  ion 


implementation  utilizing  the  VGM  LCD  is  shown  schematically  in 
Fig.  2  by  specifying  sequentially  the  nature  of  the  transformations 
from  input  intensity  to  spatial  frequency  within  the  VGM  LCD.  from 
spatial  frequency  to  spatial  Filter  amplitude  transmittance  in  the 
Fourier  transform  plane,  and  from  output  amplitude  to  output 
intensity  (usually  by  means  of  square  law  detection).  The  overall 
nonlinearity  achieved  can  be  easily  compensated  for  the  functional 
dependences  of  the  separate  steps  by  adjustment  of  the  selected 
spatial  filter  transmittance  function. 

The  VGM  LCD  has  thus  far  been  utilized  to  perform  a  wide 
variety  of  parallel  nonlinear  point  transformations,  including  level 
slicing,’  10  12  binary  logic  functions  ( AND,  OR,  NOR,  etc.),10- 1 2  and 
full  binary  addition  (inputs:  two  addend  bit  planes  and  one  carry  bit 
plane:  outputs:  sum  bit  plane  and  carry  bit  plane)."  The  purpose  of 
this  paper  is  to  describe  the  physical  principles  of  operation  of  the 
variable  grating  mode  liquid  crystal  device,  identify  areas  of  strength 
and  weakness,  and  differentiate  limitations  to  current  device  per¬ 
formance  thought  to  be  fundamental  in  origin  from  those  that  are 
seemingly  technological.  Section  2  consists  of  a  more  detailed  descrip¬ 
tion  of  the  device,  its  operating  mode,  and  its  operational  properties. 
The  fundamental  origins  of  these  operational  properties  are  examined 
in  Sec.  3,  in  which  the  natural  focus  will  be  the  physical  mechanism  of 
the  variable  grating  mode  effect  in  nematic  liquid  crystal  mixtures. 
Experimental  and  theoretical  efforts  to  elucidate  the  nature  of  this 
mechanism  are  described  in  Sec.  4,  which  concludes  with  several 
important  but  as  yet  unanswered  questions. 

2.  DEVICE  DESCRIPTION  AND  OPERATIONAL  MODE 

A  number  of  important  aspects  of  device  construction  and  device 
operation  are  reviewed  in  this  section.  A  more  complete  description 


Fig.  3.  Polarization  photomicrographs  o<  tha  liquid  crystal  domain  pat¬ 
tern  in  an  electrically  activated  call.  In  (a),  the  polarizer  was  oriented  at 
90°.  and  the  analyzer  at  90°.  with  respect  to  the  grating  wave  vector.  In 
|b),  the  polarizer  was  oriented  at  90°.  and  the  analyzer  at  10°.  with 
respect  to  tha  grating  wave  vector.  Tha  unit  vector  Hq*  denotes  tha 
direction  of  quiescent  alignment,  and  k0  indicates  tha  direction  of  tha 
grating  wave  vector. 


of  these  concepts  has  appeared  previously.’-12 

The  critical  element  of  the  VGM  LCD  is  a  thin  (4  to  12  urn)  layer 
of  nematic  liquid  crystal  mixture12  that  exhibits  a  periodic  modulation 
of  the  liquid  crystal  director,  and  hence  of  the  index  ellipsoid,  under 
application  of  an  electric  field  normal  to  the  plane  of  the  layer.  By 
means  of  suitable  preferential  alignment  techniques.’ 12  the  quiescent 
state  of  the  liquid  crystal  is  homogeneous  (parallel  to  the  plane  of  the 
layer).  As  will  be  discussed  in  Sec.  4,  this  periodic  variation  of  the 
principal  axes  of  the  dielectric  tensor  gives  rise  to  a  birefnngent  phase 
grating  characterized  by  striking  and  unique  optical  properties,14 
The  grating  can  be  visualized  in  a  polarizing  microscope,  as  shown  in 
Fig.  3,  by  utilizing  the  birefringence  properties  of  the  periodic  per¬ 
turbation.  Distinct  polarizer  analyzer  combinations  give  rise  to 
remarkably  different  grating  images,  as  can  be  seen  by  comparison  of 
Figs.  3(a)  and  3(b)  (see  Sec.  4).  Furthermore,  the  grating  period  is 
observed  experimentally  to  be  related  inversely  to  the  applied  voltage 
across  the  layer.  Above  the  threshold  for  domain  formation,  there¬ 
fore.  the  spatial  frequency  of  the  grating  is  a  linear  function  of  the 
voltage  across  the  layer,  as  shown  for  a  variety  of  nematic  liquid 
crystal  mixtures  in  Fig.  4. 

This  voliage-to-spatia!  frequency  transformation  can  be  optically 
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Fig.  4.  VQM  spatial  fraquancy  as  s  function  of  sppliad  voltaga  for  various 
namatic  liquid  crystal  mixtures. 


addressed  by  means  of  a  photoconductor  placed  in  series  with  the 
liquid  crystal  layer,  as  shown  in  Fig.  5  and  described  in  the  previous 
section.  The  photoconductive  layer  employed  in  devices  constructed 
thus  far  is  comprised  of  evaporated  or  ion-beam  sputtered  zinc  sulfide 
(ZnS),  chosen  to  optimize  the  impedance  match  with  the  liquid  crystal 
layer  (p  >  I010  ft  ■  cm).  The  layer  thicknesses  employed  were  of  order 
1.5  to  5  pm.  As  shown  in  Fig.  5,  the  photoconductive  and  liquid 
crystal  layers  are  sandwiched  between  indium  tin  oxide  ( ITO)-coated 
1.2  cm  thick  glass  optical  flats.  The  liquid  crystal  layer  thickness  is 
determined  by  a  perimeter  Mylar  spacer. 

In  operation,  a  dc  bias  voltage  is  applied  between  the  indium  tin 
oxide  electrodes,  of  order  40  to  150  V.  The  input  image  to  be  spatial 
frequency  encoded  is  focused  on  the  ZnS  photoconductor,  produc¬ 
ing  image-wise  modulation  of  the  local  voltage  across  the  liquid 
crystal  layer,  thus  effecting  a  parallel  intensity-to-spatial  frequency 
conversion.  The  high  lateral  impedance  of  the  thin  film  layers  allows 
high  resolution  images  to  be  processed  with  low  pixei-to-pixel  cross¬ 
talk.  The  device  sensitivity  is  optimized  for  exposure  at  blue  and 
near-ultraviolet  wavelengths  due  to  the  peak  photosensitivity  of  zinc 
sulfide  in  that  spectral  region.  Quasi-nondestructive  readout  can  be 
accomplished  at  wavelengths  beyond  the  photoconductivity  edge, 
such  as  that  of  the  He-Ne  laser  (6328  A).  Image  erasure  occurs  with 
removal  of  the  input  image,  within  the  dielectric  and  liquid  crystal 
relaxation  times  of  the  device  (see  Sec.  3).  To  date,  all  VGM  LCDs 
that  we  have  constructed  have  been  designed  for  transmissive  read¬ 
out.  although  reflective  readout  is  possible  with  incorporation  of  an 
appropriate  dielectric  mirror.  Such  a  configuration  would  have  the 
advantage  of  fully  separating  the  reading  and  writing  functions, 
allowing  for  increased  effective  optical  gain. 


Fig.  6.  Schematic  diagram  of  the  VQM  liquid  crystal  davica.  Currant 
davicea  arc  road  out  in  tranamiaaion  at  a  wavelength  of  photoconductive 
inaensitivity. 


3.  FUNDAMENTAL  ORIGINS  OF  THE  OPERATIONAL 
PROPERTIES 

As  was  mentioned  in  Sec.  2,  the  operational  properties  of  the  VGM 
LCD  are  primarily  determined  by  the  variable  grating  mode  effect 
exhibited  by  the  nematic  mixture  liquid  crystal  layer.  The  current  state 
of  knowledge  concerning  the  physical  origin  of  this  unique  effect  is 
summarized  in  Sec.  4.  In  this  section,  we  outline  a  number  of  key 
factors  and  considerations  that  affect  several  important  device  prop¬ 
erties  in  order  to  provide  both  focus  and  a  frame  of  reference  for  the 
succeeding  section.  These  device  properties  include  the  accessible  range 
of  spatial  frequencies,  the  number  of  accessible  gray  levels,  the  func¬ 
tional  dependence  of  diffraction  efficiency  on  applied  voltage,  maxi¬ 
mum  diffraction  efficiency,  response  time,  device  uniformity,  device 
input  sensitivity,  and  device  operational  lifetime. 

The  accessible  range  of  spatial  frequencies  extends  from  the  thresh¬ 
old  for  grating  formation  at  the  low  end  to  the  onset  of  dynamic 
scattering  induced  by  high  electric  fields  at  the  high  end.  as  shown  in 
Fig.  4.  For  phenyl  benzoate  mixtures  with  slightly  negative  dielectric 
anisotropy  ( <— 0.30)  such  as  HRL  2N40.15  this  range  extends  from 
approximately  200  line  pairs  /  mm  to  over  600  line  pairs;  mm.  In  order 
to  avoid  overlap  of  higher  diffracted  orders  from  lower  spatial  fre¬ 
quencies  with  lower  diffracted  orders  from  higher  spatial  frequencies, 
the  maximum  range  that  can  be  processed  (uniquely  assigned  to  spe¬ 
cific  gray  levels)  spans  a  factor  of  two  in  spatial  frequency.  For  exam¬ 
ple,  a  usable  range  in  HRL  2N40  extends  from  300  line  pairs i  mm  to 
600  line  pairs i  mm  without  order  overlap.  This  accessible  range  can  be 
extended  by  an  additional  factor  of  two  by  utilizing  the  orthogonal 
polarization  behavior  of  alternating  diffracted  orders,  as  described  in 
Sec.  4.  Hence,  the  accessible  range  of  spatial  frequencies  observed  in 
several  of  the  nematic  liquid  crystal  mixtures  tested  so  far  is  sufficient 
for  optimized  gray  scale  processing.  It  should  be  noted  that  although 
the  maximum  number  of  resolution  elements  that  can  be  processed  is 
linearly  proportional  to  the  highest  spatial  frequency  utilized  for 
devices  of  a  given  size  (see  discussion  below),  use  of  significantly  larger 
spatial  frequencies  begins  to  place  stringent  requirements  on  the 
Fourier  transform  lens  due  to  f-number  reduction.  For  example,  to 
utilize  a  spatial  frequency  range  of  600  line  pairs  mm  requires  the 
output  optics  to  have  an  f-number  less  than  1.2  (or  less  than  2.6  if  the 
lens  is  displaced  off-axis  to  accept  only  the  positive  diffracted  orders). 
This  is  primarily  a  pragmatic  limitation  rather  than  a  fundamental  one. 
as  VGM  effects  have  been  observed  at  spatial  frequencies  exceeding 
1000  line  pairs/ mm. 14 

The  number  of  accessible  gray  levels  that  lead  to  well-separated 
diffraction  orders  in  the  filter  plane  is  limited  by  the  ratio  of  the 
frequency  range  between  VGM  harmonics  to  the  object  spectrum 
bandwidth,  as  shown  in  Fig.  6.  The  object  spectrum  bandwidth  is  in 
turn  limited  primarily  by  two  effects;  spot  size  due  to  diffraction  from 
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Fig.  6  Gray  level  resolution.  The  numbar  o<  accaaaibia  gray  lavala  ia 
limit  ad  by  tha  ratio  of  the  aaparation  batwaan  VGM  harmonica  to  tha 
obiact  xpectrum  bandwidth. 


finite-sized  pixel  apertures,  and  grating  imperfections  that  cause 
local  deviations  from  uniform  spatial  frequency.  The  first  effect  is 
fundamental  and  has  been  treated  previously.10  The  principal  result 
of  this  analysis  is  the  inequality 


Fig.  7.  Polarization  photomicrograph  of  a  VGM  liquid  cryatal  layar  exhib¬ 
iting  a  number  of  grating  diacontinuitiea  (circled). 


bv0  >2N  .  (1) 

in  which  b  is  the  pixel  width.  v0  is  the  lowest  usable  VGM  spatial 
frequency,  and  N  is  the  desired  number  of  distinguishable  gray  levels. 
This  inequality  requires  that  the  pixel  size  contain  at  least  2N  periods 
of  the  lowest  grating  frequency  if  N  gray  levels  are  to  be  processed. 
For  example,  a  256X256  pixel  image  could  be  processed  with  32 
distinguishable  gray  levels  on  a  50  mm  square  device  with  v0  =  300 
line  pairs .  mm.  An  additional  restriction  arises  from  scattering  effects 
due  to  grating  imperfections,  which  tend  to  further  increase  the  size 
of  the  diffracted  orders.  The  most  common  type  of  imperfection 
observed  in  these  devices  is  the  joining  or  splitting  of  grating  lines,  as 
shown  in  the  photomicrograph  in  Fig.  7.  The  origin  of  these  “discli- 
nations"  is  not  at  present  understood,  although  the  density  of  occur¬ 
rence  of  such  imperfections  is  directly  related  to  the  quality  of 
substrate  preparation. 

A  typical  measurement  of  the  functional  dependence  of  diffrac¬ 
tion  efficiency  on  the  applied  voltage  across  the  VGM  liquid  crystal 
layer  is  shown  in  Fig.  8.  Since  the  applied  voltage  is  linearly  related  to 
the  induced  grating  spatial  frequency.,  this  relationship  is  illustrative 
of  the  dependence  of  the  diffraction  efficiency  on  spatial  frequency  as 
well.  To  first  order,  the  nature  of  this  dependence  is  not  important  to 
the  implementation  of  optical  processing  functions  since  any  varia¬ 
tion  in  diffraction  efficiency  with  spatial  frequency  can  be  linearized 
by  insertion  of  an  appropriate  multiplicative  filter  in  the  focal  plane 
of  the  Fourier  transform  lens.  In  any  case,  the  theoretical  functional 
dependence  can  be  derived  only  from  knowledge  of  the  relationship 
between  the  induced  orientational  angles  of  the  liquid  crystal  direc¬ 
tor  and  the  applied  voltage  across  the  layer.  This  relationship  is 
discussed  further  in  the  subsequent  section. 

The  maximum  diffraction  efficiency  that  can  be  achieved  at  a 
given  spatial  frequency  and  applied  voltage  depends  fundamentally 
on  the  magnitude  of  the  anisotropy  in  the  index  of  refraction  ( An  “ 
ne  —  n0.  with  ne  the  extraordinary  refractive  index  for  polarization 
parallel  to  the  molecular  axis  and  n0  the  ordinary  refractive  index  for 
polarization  perpendicular  to  the  molecular  axis),  the  magnitude  of 
the  penodic  angular  reorientation  of  the  liquid  crystal  director,  and 
the  thickness  of  the  VGM  liquid  crystal  layer.  Full  penodic  reonen- 
tation  of  the  index  ellipsoid  from  homogeneous  (parallel  to  the 
substrate)  alignment  to  homeotropic  (perpendicular)  alignment  for 
HRL  2N40(  An  =  0.1 5)  in  a  6  urn  cell  read  out  at  6328  \  gives  rise  to 
an  optical  phase  modulation  of  approximately  9  rad.  Hence,  the 
maximum  diffraction  efficiency  is  fundamentally  limited  to  that 
expected  for  a  pure  sinusoidal  phase  grating.1’  In  practice,  full 
reorientation  is  typically  not  achieved  before  the  onset  of  dynamic 
scattering,  although  reorientation  angles  of  45s  are  thought  at  pres¬ 
ent  to  be  commonly  reached  (see  Sec.  4).  As  can  be  seen  from  Fig.  8, 
typical  second-order  diffraction  efficiencies  are  of  order  20%.  This 


Fig.  8.  Diffraction  efficiency  as  a  function  of  applied  voltage  acroia  a  I 
nematic  liquid  cryatal  mixture  of  phenyl  benzoates  (HRL  2N40).  The  lever  I 

thickness  was  approximately  6  «m.  as  defined  by  a  perimeter  Mylar 
spacer. 


order  is  larger  than  the  first  diffracted  order  due  to  the  peculiar 
nature  of  the  birefnngent  phase  grating  formed  by  the  VGM  distor¬ 
tion  (details  are  given  for  this  phenomenon  in  Sec.  4) 

The  response  time  of  the  cell  is  a  critical  parameter  that  directly 
affects  the  achievable  overall  processing  throughput  rate.  At  present 
it  is  the  major  factor  inhibiting  widespread  incorporation  of  VGM 
devices  in  optical  processing  systems.  The  rise  time  for  grating  for¬ 
mation  from  below-to  above  threshold  varies  from  mixture  to  mix¬ 
ture.  but  is  typically  of  order  one  second.  The  response  time  for 
grating  change  in  response  to  a  step  increase  in  applied  voltage 
(corresponding  to  a  step  increase  in  grating  spatial  frequency)  is 
typically  of  order  a  fraction  of  a  second. The  response  time  of  the 
photoconductive  voltage  division  across  the  liquid  crystal  layer  is  not 
a  significant  factor  by  orders  of  magnitude  relative  to  the  reonen- 
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tation  response  time,  so  that  eventual  improvements  in  VGM  LCD 
response  time  will  accrue  only  by  advances  in  the  state  of  understand¬ 
ing  of  the  physical  origin  of  the  VGM  effect  and  the  dynamical  nature 
of  the  grating  reorientation  process,  followed  by  appropriate  modifi¬ 
cation  of  the  device  operational  mode  to  enhance  the  rate  of  molecu¬ 
lar  reorientation  and  /  or  a  search  for  a  nematic  liquid  crystal  mixture 
with  physical  characteristics  optimized  for  dynamic  VGM  effects. 

The  input  sensitivity  of  the  VGM  LCD,  defined  as  the  input 
(writing)  intensity  per  unit  area  per  unit  change  in  grating  spatial 
frequency,  is  determined  by  a  number  of  factors.  These  include  the 
slope  of  induced  grating  spatial  frequency  as  a  function  of  applied 
voltage  for  the  particular  nematic  liquid  crystal  mixture  employed, 
the  wavelength  dependence  of  the  photoconductive  layer  photosensi¬ 
tivity,  and  the  cell  switching  ratio  (fractional  increase  in  voltage 
across  the  liquid  crystal  layer  from  illumination  at  the  threshold  for 
grating  formation  to  saturation).  The  first  factor  (liquid  crystal 
response  slope)  varies  significantly  from  mixture  to  mixture  (see  Pig. 
4).  Of  the  VGM  nematic  liquid  crystal  mixtures  investigated  to  date, 
HRL  2N40  has  proved  to  be  nearly  optimum  in  this  regard.  It  is  not 
yet  clear  what  fundamentally  influences  and  eventually  limits  this 
parameter.  The  photosensitivity  of  the  photoconductive  layer  is 
determined  primarily  by  the  choice  of  photoconductive  material 
(limited  to  those  that  can  be  appropriately  impedance-matched  to 
the  liquid  crystal  layer),  method  and  quality  of  thin  film  deposition, 
layer  thickness,  spectral  width  and  central  wavelength  of  the  expo¬ 
sure  (writing)  illumination,  and  the  operational  bias  voltage  em¬ 
ployed.  The  cell  switching  ratio  is  a  function  of  the  senes  impedance 
of  the  liquid  crystal  layer,  the  impedance  of  the  unilluminated  photo¬ 
conductive  layer,  and  the  impedance  of  the  photoconductive  layer 
under  saturation  illumination.  In  addition,  the  cell  switching  ratio 
will  be  altered  by  incorporation  of  surfactant  layers  to  improve  liquid 
crystal  quiescent  alignment,  and  of  a  dielectric  mirror  in  the  reflective 
readout  device  structure.  At  this  stage  of  the  device  development,  the 
input  sensitivity  of  the  VGM  LCD  has  not  been  optimized.  A  typical 
value  of  15  (m W;  cm2)  •  (mm ;  line  pair)  was  obtained  with  a  VGM 
LCD  consisting  of  a  6  Mm  layer  of  HRL  2N40  in  series  with  a  5  Mm 
thick  evaporated  ZnS  layer  that  had  been  polished  and  rubbed  with 
surfactant  polyvinyl  alcohol,  operated  at  160  V  dc,  and  illuminated 
in  the  passband  410  to  550  nm.* 

The  uniformity  of  VGM  LCD  response  depends  inherently  on 
technological  issues,  including  uniformity  of  layer  thicknesses,  homo¬ 
geneous  mixing  of  the  liquid  crystal  material  employed,  and  the 
as-deposited  spatial  dependence  of  photoconductive  sensitivity. 
Whereas  it  is  relatively  straightforward  to  construct  electrically  acti¬ 
vated  VGM  cells  (see  Sec.  4)  that  exhibit  a  high  degree  of  spatial 
uniformity,  deposition  of  a  photoconductive  layer  with  equivalent 
spatial  homogeneity  has  proven  more  difficult.  Nonuniformity  of  the 
device  response  characteristic  can  be  a  contributing  factor  in  the 
establishment  of  the  maximum  number  of  accessible  gray  levels 
discussed  previously.  In  experimental  devices  constructed  thus  far, 
device  uniformity  has  not  proven  to  be  the  limiting  factor.  In  any 
case,  it  is  expected  that  response  nonuniformities  can  be  minimized 
significantly  by  improvements  in  the  photoconductive  layer  deposi¬ 
tion  process. 

The  lifetimes  of  experimentally  constructed  VGM  LCDs  have 
ranged  from  less  than  a  week  to  over  a  year.  The  causes  of  VGM 
device  failure  have  not  yet  been  extensively  studied,  although  several 
contributing  factors  can  be  identified.  These  factors  include  the 
purity  and  composition  of  the  liquid  crystal  mixture  employed,  the 
nature  of  the  liquid  crystal/ photoconductive  layer  interface,  the 
integrity  of  the  device  sealing  process,  and  the  device  operational 
history.  Since  the  VGM  effect  requires  a  dcapplied  voltage,  unidirec¬ 
tional  ion  poisoning  may  contribute  to  gradual  device  degradation. 

4.  PHYSICAL  ORIGIN  OF  THE  VARIABLE 
GRATING  MODE  EFFECT 

As  can  be  clearly  understood  from  the  discussion  presented  in  Sec.  3,  a 
vast  majority  of  the  important  device  operational  properties  depend 
critically  on  the  detailed  nature  of  the  grating  formation  and  dynamic 


Fig.  9.  Variable  grating  mod*  liquid  crystal  teal  geometry  showing  the 
Cartesian  coordinate  system  referred  to  in  the  text  as  well  as  the  molecu¬ 
lar  orientation  angles  a  and  n  This  configuration  was  utilizad  in  the 
polarized  light  diffraction  efficiency  and  photomicroecopy  experiments 


Fig.  1 0.  The  polarization  behavior  of  VGM  diffracted  orders,  with  illumi¬ 
nation  normal  to  the  plane  of  the  liquid  crystal  layer.  The  left-hand  column 
indicates  the  input  polarization  associated  with  each  row  of  output  polari¬ 
zations.  The  inset  shows  the  corresponding  orientation  of  the  VGM 
grating. 


reorientation  process.  The  elucidation  of  the  fundamental  nature  of 
the  grating  (molecular  orientation  angles  as  a  function  of  applied 
voltage  across  the  liquid  crystal  layer)  and  of  the  physical  mechanism 
that  gives  rise  to  the  observed  periodic  instability  and  its  reorienta¬ 
tion  dynamics  has  been  a  subject  of  considerable  experimental  and 
theoretical  interest.11  I6I,'3°  In  this  section,  the  present  state  of 
understanding  of  the  VGM  effect  is  described,  and  a  number  of 
important  unresolved  questions  are  presented. 

The  grating  produced  by  the  periodic  spatial  reorientation  of  the 
nematic  liquid  crystal  molecules  is  quite  unusual,  giving  rise  to 
striking  polarization-dependent  properties. 1 '  These  diffraction  effects 
can  be  investigated  in  an  electrically  activated  VGM  cell  with  no 
intervening  photoconductive  layer,  as  shown  in  Fig.  9.  The  orienta¬ 
tion  of  the  grating  is  such  that  the  grating  wave  vector  is  perpendicu¬ 
lar  to  the  direction  of  unperturbed  alignment,  which  is  homogeneous 
and  induced  by  unidirectional  rubbing  or  ion-beam  milling.  That  is, 
the  periodic  modulation  direction  is  perpendicular  to  the  initial  (zero 
applied  bias)  liquid  crystal  director  (long  molecular  axis),  as  shown 
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in  (he  polanzation  micrographs  (Fig.  3). 

For  all  linear  input  polarization  angles,  the  even  and  odd  diffrac¬ 
tion  orders  are  found  to  be  essentially  linearly  polanzed.  In  addition, 
the  even  diffraction  orders  are  nearly  linearly  polanzed  parallel  to 
the  “domains"  compnsmg  the  VGM  grating,  as  shown  in  Fig.  10.  For 
input  polanzation  perpendicular  to  the  domains,  the  even  orders  are 
found  to  be  almost  fully  extinguished.  On  the  other  hand,  the  odd 
diffraction  orders  are  nearly  linearly  polanzed  with  a  major  axis  that 
rotates  counterclockwise  at  the  same  rate  as  the  input  polanzation  is 
rotated  clockwise.  This  effect  is  the  same  as  that  produced  by  a 
half-wave  plate  oriented  at  43s  with  respect  to  the  grating  wave 
vector.  For  input  polarization  at  43s  to  the  wave  vector,  all  orders  are 
observed  in  the  far-field  diffraction  pattern.  An  analyzer  placed  on 
the  output  side  of  the  VGM  device  can  be  rotated  to  extinguish  the 
even  orders  (when  onented  parallel  to  the  grating  wave  vector)  or  the 
odd  orders  (when  onented  at  —45°  to  the  grating  wave  vector). 

These  unusual  polarization  properties  have  recently  been  utilized 
to  determine  the  spatial  distnbution  of  the  molecular  onentation 
within  the  VG  M  liquid  crystal  layer. 1 1  The  polanzation  dependence 
of  the  diffraction  phenomena  is  directly  related  to  the  formation  of  a 
birefnngent  phase  grating,1’  in  which  the  pnncipal  axis  of  the  index 
ellipsoid  vanes  penodically  both  in  the  plane  of  the  grating  (charac- 
tenzed  by  an  onentation  angle  a)  and  normal  to  the  plane  of  the 
grating  (charactenzed  by  an  onentation  angle  p).  The  angular  coor¬ 
dinates  are  as  shown  in  Fig.  9. 

The  polanzation  properties  of  light  diffracted  by  the  liquid  crystal 
birefnngent  phase  grating  can  be  summarized  by  means  of  a  transfer 
matrix  that  connects  the  output  polanzation  at  each  point  <  x.y)  on 
the  rear  surface  of  the  liquid  crystal  layer  with  the  input  polarization 
at  the  front  surface  of  the  liquid  crystal  layer.  On  the  basis  of  the 
expenmental  observations,  this  matnx  must  be  of  the  form 


A0*A(x:p)  B(x:2p) 

C(x.2p)  D0  +  D(x;p)_ 


(2) 


in  which  the  notation  A(  x;p)  indicates  that  the  complex  amplitude  A 
varies  in  the  x  direction  with  penodic  repetition  distance  p.  For 
uniaxial  liquid  crystal  molecules  at  an  arbitrary  onentation  (a.p). 
assumed  uniform  throughout  the  layer  thickness  at  a  given  coordi¬ 
nate  in  the  x  direction,  the  Jones  matnx  can  be  determined  by 
appropriate  rotations  of  the  index  ellipsoid,  which  yields 


2rrx 

a  =  o0  cos  . 


2  irx 

n  =  r  n0  sm  — — 


(4) 


This  procedure  allows  the  extraction  of  the  maximum  onentational 
excursion  angles  a0(V)  and  p0(V)  as  functions  of  the  applied  bias 
voltage  above  the  threshold  for  grating  formation,  as  shown  in  Fig. 
I2(a0)and  Fig.  13(p0).  In  each  case,  subject  to  the  assumed  forms  of 
a  and  p  implicit  in  Eq.  (4),  it  is  observed  that  the  maximum  excursion 
angles  both  in  and  out  of  the  plane  of  the  grating  seem  to  increase  as 
the  loganthm  of  the  applied  voltage. 

The  spatial  distnbution  of  the  ends  of  the  liquid  crystal  molecules 
descnbed  by  Eq.  (4)  is  approximately  cycloidal.  Such  a  dependence  of 
the  angles  a  and  p  on  the  spatial  coordinate  x  has  been  predicted  by 
direct  minimization  of  the  free  energy  in  a  similar  nematic  liquid  crystal 
system.21  This  panicular  solution  is  obtained  by  incorporation  of  the 
converse  flexoelectnc  effect  in  the  expression  for  the  free  energy. 

The  flexoelectric  effect  descnbes  a  strain-induced  polanzation  that 
anses  due  to  molecular  shape  effects  in  conjunction  with  a  nonzero 
dipole  moment,  as  shown  schematically  in  Fig.  14.  The  converse 
flexoelectnc  effect  thus  pertains  to  a  polanzation-induced  strain  within 
the  liquid  crystal  layer,  which  can  result  in  a  penodic  molecular  re¬ 
orientation  characterized  by  a  linear  dispersion  relation  between  the 
grating  wave  vector  and  the  applied  field,  as  is  observed  expenmen- 
tally.  Including  the  dielectnc.  distonion.  and  flexoelectnc  contnbu- 
tions  to  the  free  energy,  yields  an  expression  of  the  form 


VGM 


-/ 


=  F, 


^dielectric  *  ^distortion  ^flexoelectnc  j  dV 


»  -  j” /is  ~  *o)  <£  '  'n):dV 

lf[ 

■+•  Kj(n  x  7  x  firj 


K,(7  ■  n):  IC-(h  ■  7  x  n); 


1  -  sin*a<  1  -  e2*) 

smacosaO  —  el*) 

-  J  |e,(7  •  n)  (n  •  £) 

sinacosa(  1  —  el*) 

1  -  cosJa(  I  —  eJ*)_ 

(3) 

-  e3  j[f7  X  n)  X  n]  ■  £  j 

dV  .  (5) 


where 


in  which  t  is  the  liquid  crystal  layer  thickness,  n0  is  the  ordinary  index 
of  refraction,  n,  is  the  extraordinary  index  of  refraction,  and  K  is  the 
wavelength  of  readout  illumination  employed.  The  angles  a  and  p 
are  assumed  to  be  penodic  functions  of  x  and  independent  of  y  and  z. 
Measurement  of  the  intensities  in  each  diffraction  order  for  a  min¬ 
imum  set  of  polanzer  analyzer  onentations  uniquely  determines  the 
magnitudes  of  the  Founer  components  of  the  polanzation  transfer 
matnx.  These  expenmentally  denved  values  can  then  be  compared 
with  the  theoretically  calculated  coefficients  of  Eq.  (3)  (by  harmonic 
expansion)  for  different  possible  assumptions  concerning  the  spatial 
distnbution  of  the  orientation  angles  a  and  p.  An  example  of  such  a 
companson  between  theory  and  expenment  is  shown  in  Fig.  1 1 . 
under  the  assumption  that  the  spatial  dependences  of  a  and  p  are 
given  by 


in  which  K,.  K;,  and  K,  are  the  elastic  constants  for  splay,  twist,  and 
bend  deformations,  respectively;  n  is  the  liquid  crystal  director;  £  is 
the  applied  electnc  field:  te  and  e0  are  principal  components  of  the 
dielectnc  tensor  of  the  liquid  crystal:  and  e,  and  e,  are  flexoelectnc 
coefficients.’1  Minimization  of  FVCM  with  respect  io  the  onentation 
angles  a  and  rj  of  the  director,  subject  to  fully  pinned  boundary 
conditions  at  both  substrate  surfaces,  and  with  the  simplifying 
assumption  that  K ,  =  =  It.  yields21 


a  =  o0cos(kx)cos(  »rz  t)  . 
P  =  p0sin(kx)cos(  >rz  t)  . 


(6) 


in  which  k  is  the  grating  wave  vector,  and  t  is  the  liquid  crystal  layer 
thickness.  This  solution  generates  a  dispersion  relation  between  E 
and  k  of  the  form 


’k:  + 

■  (>r  t);_ 

2 

k: 

+  M 

_k*  i-  ( ir/t)*_ 
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Fig.  11.  Measured  diffracted  order  intensities  as  a  function  of  theoretical 
intanaitiaa  caiculatad  from  tha  uniaxial  VQM  modal  daacribad  in  tha  taxt. 


Fig.  13.  Tha  out-of-plana  molecular  orientation  angle  (n0)  aa  a  function  of 
tha  applied  dc  biaa  voltage  arroaa  tha  call  (V)  (aee  Fig.  5). 


Fig.  12.  Thein-planemolecularorientationanglelao)  at  a  function  of  tha 
applied  dc  biaa  volt  age  acroea  tha  call  (V)  (aaa  Fig.  5). 


in  which  e*  ™  e,  —  e3,  u  “  (»aK.  4Tre*2),  and  ta  “«e  _  f0.  The 
dispersion  relation  is  linear  when  k  »  rr  t,  as  shown  in  Fig.  15 
(compare  with  the  expenmental  relationship  shown  in  Fig.  4). 

The  origin  of  the  periodic  instability  in  VGM  liquid  crystals  has 


occurrence  of  tha  flexoelectric  affect  due  to  a  shape  anisotropy  in  liquid 
crystal  molecules  with  a  permanent  dipole  moment.  In  (a),  the  molecular 
orientation*  are  randomly  distributed,  resulting  in  taro  net  polarization.  In 
(b),  the  applied  distortion  induces  a  shape-dependent  molecular  realign¬ 
ment.  resulting  in  a  net  polarization.  (After  Ref.  31  I 


not  yet  been  established  beyond  doubt,  although  the  accumulated 
evidence  points  strongly  toward  the  converse  flexoelectric  effect. 
This  assignment  is  also  intuitively  appealing  since  the  VGM  effect  is 
observed  only  in  liquid  crystal  mixtures  with  slightly  negative  dielec¬ 
tric  anisotropy,12  for  which  both  the  dielectric  and  distortion  contri¬ 
butions  to  the  free  energy  increase  for  deviations  from  uniform 
alignment.  The  addition  of  the  flexoelectric  term  counteracts  these 
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mm  COM— Bonding  to  HRL  2N40  pbwtyl  benzoate  mixture  [aoa  Eg.  (7)1. 
subject  to  the  aaaumption  that  K,  =  K2  =  K. 


effects,  producing  a  free-energy  minimum  at  nonzero  distortion 
angles.  Furthermore,  this  model  predicts  a  static  periodic  perturba¬ 
tion:  no  hydrodynamic  or  electrohydrodynamic  effects  are  observed 
in  our  test  cells. 

The  dynamics  of  grating  formation,  reorientation,  and  relaxation 
are  subjects  of  ongoing  experimental  and  theoretical  investigation. 
Understanding  of  the  basic  physical  principles  underlying  these 
effects  is  vital  to  the  success  of  efforts  to  improve  the  VGM  LCD 
response  time. 

A  number  of  spatial  discontinuities  which  are  similar  in  appear¬ 
ance  to  crystallographic  dislocations  can  be  observed  within  the 
liquid  crystal  layer  (see  Fig.  7).  The  "dislocation"  density  increases 
with  time  as  each  VG  M  device  deteriorates,  and  strongly  depends  on 
the  manner  in  which  the  applied  bias  is  brought  from  below  to  above 
the  threshold  of  grating  formation.  In  some  cases,  as  the  applied  bias 
across  an  electrically  activated  cell  is  increased  or  decreased,  the 
dislocations  propagate  past  each  other  while  the  grating  period 
decreases  or  increases,  respectively  Dislocations  with  opposite  orien¬ 
tations  propagate  in  opposite  directions  until  a  new  stable  equilib¬ 
rium  situation  is  achieved.  It  is  not  yet  clear  whether  these  disloca¬ 
tions  represent  true  disclinations.  for  which  the  liquid  crystal  director 
is  discontinuous  near  the  defect,  or  whether  they  are  in  fact  alterna¬ 
tive  local  continuous  solutions  to  the  free-energy  minimization,  per¬ 
haps  induced  by  an  as  yet  undiscovered  perturbation. 

A  large  number  of  research  directions  have  been  described  in  the 
foregoing  paragraphs.  In  addition,  several  others  are  worthy  of  note. 
First,  numerous  experimental  measurements  of  the  off-diagonal 
elements  in  the  Jones  matrix  describing  polarized  light  propagation 
through  the  VGM  cell  have  revealed  a  significant  asymmetry  not 
predicted  by  the  uniaxial  model.  The  origin  of  this  “B  C”  asymmetry 
effect  has  not  yet  been  elucidated.  Second,  the  thickness  dependence  of 
the  molecular  orientation  angles  a  and  n  has  not  been  fully  established. 
The  solution  proposed11  in  Eq.  (6)  assumes  full  surface  pinning  at  the 
substrate  boundaries,  and  represents  the  lowest  order  z-dependem 
mode.  The  polarization-dependent  diffracted  order  measurements 
described  herein  do  not  provide  clear  differentiation  between  a  uni¬ 
form  z-dependence  and  the  lowest  order  mode.  Third,  the  nature  of 
the  transient  solutions  to  the  free-energy  minimization  incorporating 
molecular  dynamic  and  viscosity  effects  has  not  been  treated.  Solu¬ 
tion  of  this  problem  is  key  to  response  time  optimization  of  the  VGM 
liquid  crystal  device.  Fourth,  it  is  not  yet  clear  whether  a  nematic 
liquid  crystal  mixture  can  be  found  that  exhibits  a  useful  VGM  effect 
under  ac  bias.  Such  a  mixture  may  provide  significant  immunity 
from  long-term  ion-poisoning  effects.  Finally,  it  should  be  empha¬ 
sized  that  the  variable  grating  mode  liquid  crystal  effect  provides 


only  one  possible  means  of  achieving  parallel  intensity-to-position 
encoding.  The  processing  potential  of  the  intensitv-to-position  algo¬ 
rithm  should  provide  more  than  adequate  inducement  to  intensify 
the  search  for  alternative  implementations 
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